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Abstract 



New exact upper and lower bounds are derived on the spectrum of the square of the 
hermitian Wilson Dirac operator. It is hoped that the derivations and the results will be 
of help in the search for ways to reduce the cost of simulations using the overlap Dirac 
operator. The bounds also apply to the Wilson Dirac operator in odd dimensions and are 
therefore relevant to domain wall fermions as well. 
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Introduction 

Let D{m) denote the continuum Euclidean Dirac operator where the real parameter m 
is the fermion mass. In even dimensions d a generalization of 75 exists and shall be denoted 
by 7d+i. D{0) is antihermitian and anticommutes with 7^+1. Then, H{m) = 7d+i-D(m) is 
hermitian. H{m) will be referred to as the hermitian Dirac operator. A characteristic 
property of this operator is the range of its spectrum as a function of the real mass 
parameter m. Since H^{m) = D"^ [m)D{m) it is meaningful to consider the spectrum 
of H'^{m) both in even and odd dimensions. 

Figure 1 displays the familiar spectral structure of H{m) in the continuum in an 
arbitrary fixed gauge background. The boundaries shown come from rigorous lower bounds 
on the spectrum of H'^{m). These bounds hold for any gauge background and are often 
saturated, for example in the case that the gauge background is trivial, or in the case that 
it consists of a gauge field carrying non-zero topology. There is no upper bound on i7^(m), 
and the spectrum will indeed increase indefinitely in any fixed smooth gauge background. 
All this holds also on a compact manifold, henceforth taken to be a flat torus. 
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Figure 1 Spectrum of the Dirac Hamiltonian in the continuum. All oblique lines have 
slopes ±1. 

The objective of this paper is to clarify what happens when the massive Dirac Hamil- 
tonian is put on the lattice following Wilson's prescription. The most fundamental feature 
of a lattice operator is that its spectrum is absolutely bounded from above - this is how 
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the lattice acts as a regulator. However, lower bounds obeyed by the hermitian Wilson 
Dirac operator, i7^(m), are also very important, because often we wish to use HT^{m) to 
put massless, or almost massless quarks on the lattice. 

When the gauge background is trivial, Hwijn) can be explicitly diagonalized and one 
finds the spectral structure shown in Figure 2. A simpler derivation is contained in what 
follows. 
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Figure 2 Spectrum of the Wilson Dirac Hamiltonian on the lattice for d=4. All oblique 
lines have slopes ±1. 

When the gauge field is turned on the figure gets distorted. The upper bound on 
H^{m) remains unchanged, and so does the lower bound for positive values of the mass 
parameter m. Changes occur only for m < and for the lower bound of H^{m). So long 
we are close to the trivial case the distortion is small: it amounts to the replacement of 
the string of rhombi in Figure 2 by a string of smaller rhombi, inscribed into the ones we 
have in Figure 2. The new rhombi no longer touch each other. As m is varied, eigenvalues 
of Hy^rijn) can cross zero in the intervals that open up, separating the rhombi. When the 
gauge background is random enough the internal rhombi close up completely and very low 
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eigenvalues of H^(m) are no longer excluded [1] for any mass in the segment (— 2d, 0). 
For any gauge background the figure stays mirror symmetric about the m = —d vertical 
line. 

Although we focus on even dimensions here, so long we phrase the results for the 
Wilson Dirac operator itself and not its hermitian version, they hold for odd dimensions as 
well. In particular, the five dimensional case applies to domain wall formulations of QCD. 

Notations and Conventions 

Let us start by establishing our notation. We are working on a d-dimensional hyper- 
cubic lattice. When comparing to the continuum the lattice spacing is denoted by a. On 
its links we have SU{n) matrices U^{x) which make up the gauge background the fermions 
interact with, ji = 1, 2, ... d denotes positive directions and x denotes a lattice site. The 
lattice is finite. 

The fermions are vectors V'a(^)- o; is a spinorial index, i is a gauge group index and 
X is a lattice site. The action on the fermions is described in terms of several unitary 
operators. First are the Euclidean Dirac 7^'s which act only on spinorial indices. Second 
come the directional parallel transporters which act on the site index and the group 
index. They are defined by: 

A third class of unitary operators implements gauge transformations, each characterized 
by a collection of g{x) e SU{n) acting on -0 pointwise, and only on the group indices. The 
action is represented by a unitary operator G{g) with {G{g)il)){x) = g{x)il>{x). The 
operators are "gauge covariant", 

G{g)T,{U)G^{g)^T,{U<^), 

where, 

U^^{x)^g{x)U^{x)g\x + il). 

The variables are distributed according to a probability density that is invariant 

under U ^ for any g. 

The lattice replacement of the massive continuum Dirac operator, D{m), is an element 
in the algebra generated by T^, T^, 7^. Thus, D{m) is gauge covariant. For U^{x) = 1 
the become commuting shift operators. 
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The Wilson Dirac operator, Dwi'm) is the sparsest possible analogue of the con- 
tinuum massive Dirac operator which obeys hypercubic symmetry. Fixing the so called 
r-parameter to its preferred value (r = 1), Dw{m) can be written as: 

D^ = m + J2{i-V,y, ViV, = l; V, = ^-^T, + ^-^T^ 

In even d we associate to the Wilson Dirac operator the hermitian Wilson Dirac operator, 
Hw{m) = 'jd+iDwim). 

All our lattices are assumed finite an therefore all our operators are finite dimensional 
matrices. An eigenvalue of a matrix A will be denoted by X{A); if the eigenvalues are 
labeled, the label is attached to A. When it makes sense, we may deal with the maxi- 
mal(minimal) eigenvalues of A, Ajnax(min)(^)- We choose the following norm definition for 
matrices A: \\A\\ = [Amax(^^^)]^- This is a standard choice, induced by the vector norm 
||f IP = where / is a generic component index [2]. The norm of a gauge covariant 

matrix is gauge invariant. 

Formal Continuum Limit 

The connection to the continuum is as follows: Assume to be given smooth functions-^! 
A^{x) on the torus. Then, 

= Pexp[i J dxnAij,{x)] (the symbol P denotes "path ordering"). 

Consider a smooth function tjjc{x) with same index structure as the corresponding object 
on the lattice. By looking at x's coinciding with a lattice point one gets a lattice vector 
iIj{x = fia), where n e Z*^. The action of the produces another lattice vector ij/. One 
can define a continuum operator, T^c such that the lattice restriction of T^c'^c will be a 
function ij)'^ whose lattice restriction is ij)' . The formula is-'^^: 

T^e = e"^^ D^^d^ + iA^. 

f'^ In general the A^(a;) aren't smooth functions, rather they make up a one form 

Aij^{x)dxi^ which is a smooth connection on a possibly nontrivial bundle with structure 

group SU{n) over the four-torus. 

This generalizes an observation of van Baal [3] . 



U,,ix) = lim 
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The simplicity of this expression can be viewed as a motivation to introduce the T^'s as 
central objects on the lattice in the first place. 
The formula is easy to prove: 

for any vector ip^- On the other hand, for any operator Oc acting pointwise by Oc{x) we 
have: 

Oc{x + hp) = e^^^Oc{x)e-^^^. 

Inserting this expression (with b = k-^) repeatedly into the definition of U^{x)^ imple- 
menting the shift of the argument of il}c{x) as above, and taking N to infinity at the end, 
produces the desired result using Trotter's formula [4] . 

The V]u's have associated continuum operators V^ci given by: 

V^c = e~^"^^^'^ no sum on fj,. 

The Wilson Dirac operator is a lattice restriction of the continuum operator 

aDwcim) = m + ^(1 - e'"'^^^'^). 

DwcinT) could be viewed as an approximation to ^^D^ in the continuum which is good for 
eigenvalues small in absolute value but whose spectrum is restricted to a bounded domain. 
Such operators are frequently introduced when one regulates infinities in the continuum. 
The continuum Dirac operator 1^J.D^ formally emerges as a goes to zero, and the mass 
is of order ^, where m is a pure number. But, as an operator in the continuum, Diyc(m) 
is special: when it acts on ipc to produce V'c) the values of ip'^. at lattice points are solely 
determined by values of '0c at lattice points. Therefore, there exists an exact relation to 
the lattice operator Dw{fn)- 

There is no remnant of chiral symmetry (for even dimension d) because DvKc("^) isn't 
just a function of Yl/nliJ-^tJ-'i only in the small a limit (strictly speaking, one would need 
to replace m by rrica before taking a to zero) do we get an expression involving only the 
chiral combination YIuxIij-^ij.- 

It is important to appreciate that one does not need Dvkc(O) to anticommute with 
^d+i to have some amount of lattice chirality: any reasonable Dc{m) that is a function of 
only the combination YlnxltJ-^iJ' would do. For example, if aDwcii^) were replaced by 

aD'^^{m) = m + 1 - e^M '^^^^ , 
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we would have enough symmetry because 



Si tier dete 



is unity ^logdet 



is odd in ji and this is 



enough to eliminate additive quark mass renormalization. However, the operator e 
cannot be restricted to the lattice because when it acts on il) and produces ip' it is not true 
that the values of if)' at lattice points depend only on values of if) at lattice points. 

One can try to "improve" Dwc{fT^) by looking at the difference D'y^^{m) — Dwcif^) 
to leading order in a and replacing it by a function of the T^c (again to leading order in 
a). Adding the new term to i5vFc(^) produces an operator which can be restricted to the 
lattice and is "clover improved"; it agrees with D'yy^^m) to leading and subleading order 
in a. In fluctuating gauge field backgrounds one changes the coefficient of the new term 
to a number determined numerically. 

One can also maintain chiral symmetry on the lattice exactly [5] [6] , using the overlap 
Dirac operator. 



Upper bound 



Our first objective is to find a bound for the largest eigenvalue of H^. Clearly, 
Amax(-f^H') = The triangle inequality then gives: 

< |m + d| +^ IIF^II = \m + d\ + d. 

The lowest upper bound as a function of mass is obtained at m = —d, which is a symmetry 
point for if^(m), because Dw{—d) and —Dw{—d) are unitarily equivalent. This is a 
consequence of the existence of a unitary hermitian operator S such that SV^S = — V^, 
implying SDw{'m)S = — Z)vk(— m — 2(i); S is diagonal and the diagonal entries are 1 if the 
site X has even and —1 otherwise. 5" exists because the hypercubic lattice we are 

working on is bipartite. 

For m > —d the upper bound is attained iff there exits a vector i(; which is a common 
eigenvector to all d operators, with the eigenvalue —1 in each case. It is likely to find 
such an eigenvector when [T^, T^] = for all /x and v. These commutators vanish when all 
plaquette parallel transporters are unity; this is so, in particular, in the free case. 



Lower bound 
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Let us introduce some shorthand notation: 



The unitarity of holds because of the identities 

hl-al = l, [hf„a^]=0. 

Let A(m) = X{Hw{iTT')) be some eigenvalue of Hwifn). A(m) is differentiable because 
Hw{rn) depends smoothly on m: = ^:r,i,a,(3i'a* i^h^ a,(i%{x), where Hw{m)ilj = 
X{m)ip and i/j has unit norm. Since 7I = 1 one has 



dX 



dm 



< 1. 



The theoretical usefulness of expressions for ^ has been recently emphasized by Kerler 
[7] . This inequality restricts the slope of lines describing the flow of eigenvalues of Hw (m) 
as a function of m. We shall refer to this inequality as the "flow inequality". It has an 
important consequence that we shall prove below: If we know that < Xinin{H^{m)) for 
some m, we have 

— \m — m \. 

Before describing the proof let us note that the result is useful only if 

\m — m!\ < [Amin(iy^(m))] ^ . 

1 

The main observation is that a lower bound on [Amin(iy^(m))] ^ at an arbitrary mass 

1 

point m can be extended to a lower bound on [Amin(i?^(m'))] ^ in some mass mass range 
around m. 

The basic inequality can be best proven appealing to a sketch shown in Figure 3: The 
graphical meaning of the inequality is that (m) has no eigenvalues in the area bounded 
by the right angle rhombus in the figure when it is given that there are no eigenvalues 
along its main diagonal {A^ B). Recognizing this, the proof becomes trivial: if we did have 
an eigenvalue anywhere inside the rhombus the flow inequality would have to be violated 
somewhere in order to avoid an eigenvalue flow crossing the main diagonal. 
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Figure 3 A rhombus containing no eigenvalues. All oblique lines have slopes ±1. 



We start with an explicit formula for H^{m): 



m 



m 



+ 2(m + 1) ^(1 -h^) + J2 - - - a^a^- - ^-]] • 



While in the continuum D^^ (m)D{m) commutes with 7^+1 the last term in H^{m) does 
not. All terms are individually hermitian. Since Hw{m) connected sites x, x' with'^^ 



\x — X 



0, 1 we could have expected H'^{m) to connect sites with \x 



X 



0,1, A 2, 



but because of the relations /i^ — = 1 and [h^,a^] = sites with |x — x'| = 2 are 
still disconnected. Another special property of H^{m) is that the site diagonal piece is 
proportional to the identity matrix. 

If [T/^, T„] = we have Zl^/i. ^z^^;^ = 0^ [^/x, a„] = and [/i^, h^] = 0. Then, 

Hl^im) =m^ + 2(m + 1) ^(1 - h^) + ^(1 - h^){l - K). 

If we keep all /i^ fixed but one, say h^, the dependence on the latter is linear, so the 
extremal values are obtained at h^, = ±1. The argument is applied again and again to 



■^3 For two sites x and y we define \x — y\ = \/J2ni^fj- ~ Vfj-) 
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a decreasing number of remaining directions leading to the conclusion that in order to 
find the extrema of H^, viewed as a function of the quantities /i^ (more precisely, their 
eigenvalues, since the /i^ can be simultaneously diagonalized by assumption) we only need 
to check the 2*^ possibilities = ±1. The upper bound comes out as above, and the lower 
bound on [Ainin(-ff^)] ^ has the shape shown in the Figure 2. We learn that at the points 
m = 0, —2, —4, . . . — 2d the theory has massless fermions; the multiplicities are given by 
(d-n)!n! '^^cre 171 = —2n, 77. = 0, 1, 2, . . . , d. Thus, for m = —2, —4, ... — 2d + 2 we have 
several doublers, the number of different species being given by the number of different /i^ 
configurations producing a zero at the respective special mass point. 

From now on we shall concentrate on the region — 2 < m < 0. This region is interesting 
when we want to deal with one Dirac fermion and avoid doublers. The region close to 
m = is important for traditional numerical QCD with Wilson fermions. The region close 
to m = — 1 is important for applications of the overlap Dirac operator where one would 
like H^{m) to have a large gap around zero. When [T^,Ti^] = 0, the highest lower bound 
is obtained at m = — 1. As long as all operators [T^,Ti^] are small in norm we expect the 
same to be true. We therefore focus on the point m = — 1 first, and later extend the bound 
to a range around m = — 1 using the consequence of the flow inequality established earlier. 

In the general case where the matrices do not commute, we have 

H^{-1) = 1 + J] [(1 - h^){l - K) - a^a^ - [a^, K]] . 

We now analyze each term in the bracket individually; we treat them separately because 
their spinorial index structures are different. The first term is rewritten as: 

[(1 - M(l - = i E(l - ^m)(1 - - - Tl) = Q + X 
Here, Q is positive semidefinite, 

while X depends only on T^-commutators: 

Proceeding, we find 

- E " ~^ E ^v^^^^p^ -TI^T^- TI] = Y, 
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and 

- J2 h-] = lYl - Tu] + h.c.) + (7^ + 7,)([T^, Tl] + h.c.)] = Z. 

The operators Q, X, Y, Z are all hermitian. Moreover, each of the traces of X'^, Y'^, Z'^ are 
linearly related to the single plaquette Wilson action (see below) and decrease when the 
latter increases and the continuum limit is approached. 

Consider now the commutators [T^,Tj^]. Their norm is determined by: 

[T^, T,]t[r^, T,\ = (1 - P^,)t(l - P^,), 
where the unitary P^^^y are given by: 

p _ rp'^rp'^rp rp 

The operators P^^ are site diagonal, with entries that are parallel transporters round 
plaquettes: 

{P^ytp){x) = U^jy{x)'tp{x), U^u{x) = Ul{x - i')Ul{x - u - fi)Ujy{x - u - fi)U^{x - fi). 

Unuix) is associated with the elementary loop starting at site x, going first in the negative 
V direction, then in the negative direction, and coming back round the plaquette. 
The main relation is: 

||[T^,T.]|| = ||1-P^.||. 

Any pure gauge action with the right continuum limit will strongly prefer configurations 
where all U^y{x) are close to unit matrix. Therefore, it is not unreasonable to impose the 
constraint, for all i-i > ly, 

Note that this is equivalent to 

||1 - U^y{x)\\ < 

for every site x. It is easy to see that the same bound will hold when we interchange in the 
commutator the |U, v indices, and when we replace, independently, the and operators 
by their hermitian conjugates. 

Using the triangle inequality and that ||^-B|| < ||A||||P||, we now obtain: 

/u>i/ ij,>i' jj.yi' 
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The -\/2 factor comes in because (7^ ± 7j/)^ = 2 for 7^ i/. We finally obtain: 

/lt>i^ /u>i^ p>i^ 

By the variational principle and the positivity of Q we arrive at: 

Ami„(iy'(-1)) > 1 - (2 + V2) ^ e^,. 

Our result is meaningful only when the number on the right hand side in the above equation 
is non-negative. 

In the rotational invariant case one could set e^j^ = ry. Then, for d = 4, we obtain 

VAmin(iy2(-i)) > ^1-6(2 + ^2)77 ^ yr^2o:5^. 

The general bound we obtained is: 
>^m\n{D\^{rn)Dw{rn)) ' > 
This bound is useful only for 



l-(2 + V2)^e^, 



- |1 + m|. 



11 + ml < 



1 - (2 + V2) 5^ e^. 



This range is contained in the open segment —2 < m < 0. The bound holds in both even 
and odd dimensions. In the particular case of domain wall fermions, plaquettes parallel to 
the extra dimension make no contribution since their e^i, vanishes. 



Comparison to other work 



Related issues were studied in [8] and in [9]. The authors of [8] established the upper 
bound 

in four dimensions with the restriction — 2 < m < 0. This is compatible, but less stringent 
than our upper bound, which becomes m + 8 in this mass range. 

In numerical investigations with pure gauge Wilson action, it was reported in [8] that, 
for 13 = 6.0, 6.2, 6.4 and m = -1.0, -1.2, -1.4, -1.6, for SU{3), Xma^{H^{m)) stays around 
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41 and hardly changes. Our upper bound for m = —1.6 is 6.4^ = 40.96 and increases for 
the fower m's. Thus, at the extremal mass value (assuming the value quoted in [8] was 
rounded), our bound is saturated to numerical accuracy. The claimed mass independence 
seems surprising, and not entirely consistent with numerical results at other /? values and 
volume sizes 

In [8] a bound on Xmin{H^{—l)) is also established. It is expressed in terms of a 
bound on the norm of the commutators, but the precise definition of the norm used is not 
given. I shall assume it is the one adopted in this paper. A bound is quoted only for the 
d = A case, for m = — l'^^ and for the rotational invariant case e^i, = rj. The bound derived 
in [8] is: 

Amin(^^(-l))>l- 3077. 

This bound is compatible with the result of this paper, Amin(-H^^(— 1)) > [1 — 6(2 + \/2)rf\, 
but weaker. To be sure that i7^(— 1) has no zero eigenvalues the bound in [8] places a 
restriction on -q that is stronger than ours by about one third. 

Lessons 

Let us first identify what about our results could have been expected without any 
calculations. Clearly, we know that there will be some uniform upper bounds on the 
spectrum just by virtue of compactifying momentum space and because the U^{xys are 
unitary. Moreover, once the free case is worked out and the spectral restrictions of Figure 2 
are derived, one knows that close to the continuum the structure will be essentially similar 
even in the presence of nontrivial gauge fields. The reasoning is as follows: We are dealing 
with operators that are analytic in and Figure 2 holds whenever all T^'s commute. All 
that enters in the bound derivations above is that the T^'s are unitary. Commuting unitary 
T^'s can be smoothly deformed into non-commuting ones and the changes in the spectrum 
must be smooth too. Thus, if the commutators of the T^'s are sufficiently small there will 
be a region around m = — 1 where the spectrum of H-\y{m) will have a gap around zero. 
One can simply think about the commuting case as a "semiclassical" approximation to the 
non-commuting case. 

The numerical work was carried out by the SCRI group, at the time consisting of 

Edwards, Heller and Narayanan. 

According to David Adams [9] , in unpublished work, the authors of [8] have extended 

their bound by using the triangle inequality to a range of m values contained within the 

segment (—2, 0). 
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The operators connect only sites one spacing apart in the //-direction. The gauge 
invariant norm of the commutators cannot depend on anything else but the norm of the 
elementary plaquettes. Forcing all unitary plaquette operators close to identity produces 
a link configuration for which the T^'s almost commute. The precise relation between 
the T^-commutators and the plaquettes is well known [10] since the discovery of large n 
reduction of lattice gauge theories [11]. 

So, all that really required some work was to turn the above into a quantitative 
estimate. Because of the practical difficulties associated with low eigenvalues of H^{m) 
it makes sense to try to be as careful as possible in deriving the quantitative form of the 
bounds. Still, it is known that the lower bounds in the — 2 < m < region are not directly 
useful in backgrounds generated at coupling constants that are practical in numerical QCD 
today. In spite of this, the exact bounds and their derivation might provide helpful insights, 
in particular in the context of implementations of the overlap Dirac operator. In this case 
one wishes to work with operators H^{m) with —2 < m < but with as large a gap 
around zero as possible. This would make the matrix Hwifn) well conditioned and speed 
up the calculations. 

The most basic observation is that one can control the gap in Hw{fn) by controlling 
the plaquette variables alone-'^^. This was understood long ago [13]; a natural guess would 
be that replacing the pure Wilson gauge action by the so called "positive plaquette" model 
[14] (for gauge group SU{2)) will create a gap around zero. Numerical checks by Urs Heller 
in early 1998 have shown that this was not the case [15]. In addition, one cannot just 
change the form of a single plaquette action and get something useful in four dimensions. 
The correlation length increases exponentially as the plaquettes are forced to identity 
and physical realistic volumes rapidly become totally impractical. A milder approach is 
therefore called for. There are a few possibilities. 

First, one could use a more complicated action then a single plaquette one. The idea 
is that a more complicated action might make the plaquettes close to unity, but still keep 
the gauge fields sufficiently random so that the correlation length does not exceed a few 
lattice spacings. The improvements observed in simulations using domain walls (which can 
be viewed as a particular truncation of the overlap [16]) when one switches from Wilson to 
so called "Iwasaki actions" might be a reflection of this mechanism [17]. A more systematic 
approach would be to follow an approximate renormalization group trajectory [18], where 
the correlation length is controlled, to regimes in the coupling constant space where the 

This was exploited when the parameters of the first dynamical simulation of the 
exactly massless Schwinger model were chosen [12]. 
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single plaquettes are closer to unity. A note of caution: the inclusion of the fermionic 
determinant in the gauge measure may be important and a fix that works for quenched 
simulations may fail in the dynamical case [19]. 

Another observation is that making only the plaquettes in some directions close to 
unity would help. This only requires to increase one dimension of the lattice and there is no 
exponential relation between this dimension and the closeness of the time-like plaquettes 
to unity. In four dimensions there are other good reasons for working on asymmetric 
lattices [20], so this looks like a cheap and attractive alternative worth exploring -^^ . In 
lower dimensions than four the impact of going to asymmetric lattices would be even more 
pronounced. 

Yet another possibility is to filter out the "roughness" from the gauge background 
seen by the fermions by replacing the link variables U^{x) by new link variables U^^^{x) 
which are functions of the original link variables, transform the same way under gauge 
transformations, but produce plaquette variables closer to unity. Recent work has obtained 
such "APE smeared" U^^^{x) [21] with associated plaquettes extremely close to unity 
[22]. Of course too much "filtering" may take the lattice theory at typical simulation 
parameters too far away from the desired continuum limit of QCD-^*. If this is true, one 
could also try a "half smeared" approach where only the links entering the "Wilson mass 
term" ^ J^i^C^fj- + ^/I) are APE smeared but the links entering the chiral part 

I 7^(T^— T^) are not, so the fermions are not insulated from the ultraviolet fluctuations 
in the gauge field. Unfortunately this would spoil the relations h^—a^ = 1 and [/i^, a^] = 0, 
so the consequences on the bounds are complicated. Also, the spinorial structure no longer 
only involves the projectors |(1±7^) which causes some numerical overhead. Note however 
that with APE smearing the difference between U^^^{x) and Un{x) goes to zero when the 
original commutators go to zero. Therefore, some bounds of similar structure to the 
bounds presented here would still hold. 

It is hoped that the analysis of this paper would prove helpful in guiding our search 
for improvements in the gauge action and in the structure of DwinT)- 
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